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ABSTRACT 

We report spatial fluctuation analysis of the sky brightness in near-infrared 
from observations toward the north ecliptic pole (NEP) by the AKARI at 2.4 
and 3.2 pm. As a follow up study of our previous work on the Monitor held of 
AKARI, we used NEP deep survey data, which covered a circular area of about 
0.4 square degrees, in order to extend huctuation analysis at angular scales up 
to 1000". We found residual huctuation over the estimated shot noise at larger 
angles than the angular scale of the Monitor held. The excess huctuation of the 
NEP deep held smoothly connects with that of the Monitor held at angular scales 
with a few hundreds arcseconds and extends without any signihcant variation 
to larger angular scales up to 1000". By comparing excess huctuations at two 
wavelengths, we conhrm a blue spectrum feature similar to the result of the 
Monitor held. We hnd that the result of this study is consistent with Spitzer 
Space Telescope observations at 3.6 pm. The origin of the excess huctuation 
in the near-infrared background remains to be answered, but we could exclude 
zodiacal light, dihuse Galactic light, and unresolved faint galaxies at low-redshift 
based on the comparison with mid- and far-infrared brightness, ground based 
near-infrared images. 

Subject headings: galaxies: clusters: intracluster medium - galaxies: high-redshift 
- infrared: diffuse background - methods: data analysis - stars: population III - 


zodiacal dust 
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1. INTRODUCTION 

It has been known that there exists excess extragalactic background light (EBL) 


board Cosmic Backqround Explorer (COBE 

Hauser et al. 

1 

998 


Dwek & Arendt 

1998: 

Wright & Reese 

2000; 

Cambresv et al. 

2001: 

Arendt & Dwek 

2003 

: Levenson et al 

20071 


and the Near-infrared S pectrometer (NIRS) on board Infrared Telescope in Space (IRT S; 


Matsumoto et ah 


2005h . Recent analysis of data obtained by AKARI (ITsumura et ah 


2 OI 3 I) also shows a consistent result with them. The excess EBL is several times brighter 


than the integrated light of galaxies and shows a blue stellar-like spectrum. Since the 
near-infrared sky brightness is dominated by the zodiacal light, which is scattered sunlight 
by the interplanetary dust (IPD), the absolute level of the excess EBL has uncertainty 
due to the inaccuracies in the model of zodiacal light. On the other hand, the fluctuation 
of near-infrared sky brightness could provide complemen tary information on the origin of 


ex cess emission sinc e the zodiac al lig ht is fairly smooth I Pvo 


of 


Kashlinskv et al 


(12005 


2007) and 


Matsumoto et al. 


et al. 


(12011 


2 OI 2 I ). Recent studies 


, hereafter Paper 1) using 


the data taken by space telescopes found excess fluctuations of the sky brightness at 


angular scales of a few hund reds arcseconds that c anno 


components. More recently. 


Coorav et al. 


(l2012bh and 


be explained wit 


Kashlinskv et al. 


1 known emission 


2012) found an 


almost flat angular spectrum of the fluctuation extends up to degree angular scales by using 
data set from Spitzer at 3.6 and 4.5 /rm. Characteristic features of observed fluctuation 
are a bl ue stellar-like s pectr um and a good spatial correlation between different wavelength 


bands. 


Zemcov et al. 


(120141) measured significant excess fluctuations at 1.1 and 1.6 /rm 


using the data from Cosmic Infrared Background Experiment (CIBER). 


One of the intriguing interp retations regarding the 


EBL is the Population III stars (IKashlinskv et ah 


2005 


origin o 


the near-infrared excess 
Paper il l that cannot be identihed 
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individually because of their faintness. However, there are so me questions and debates 


about such interpretation. According to 


Madau fc SilkI ( 2005 1. large excess of the EBL 


at near-infrared wavelength bands requires a large number of Population III stars in the 
early universe and produces too much heavy elements at high redshifts. Also TeV 7 -ray 


observations prefer lower ne ar-infrared EBL t 


in 7 -ray spectra of blazars (Aharonian et ah 


H.E.S.S. Collaboration 


ran observed one to explain 


2006; 


Gilmore et ah 


2012 


absorption 


Mever et al 


eatures 


2012 


201311 . In addition, recent simulations predic ted smaller fluctuatio ns 


than observational measurements by ab o ut an order of magnitu de flCoorav et al 


Fernandez et al. 


2012 


Yue et al. 


2018all 


Coorav et al. 


2012 a; 


fj 2012 bll have argued that the 


intrahalo light (IHL) due to the stars that escaped from galaxies can explain the observed 
excess fluctuation at large angular scales. 

Clearly, more observational studies are necessary to understand the nature of the 


infrared background radiation. In this stuc 


y, we extend the scope of our previous analysis 


of the AKARI Monitor held data flPaper ll) for huctuation of the sky brightness to larger 
angular scale by using the north ecliptic pole (NEP) deep survey data taken by AKARI 


(iMatsuhara et al. 


200611 . Compared to the previous study with the Monitor held, the NEP 


deep held has shorter exposure time but larger areal coverage, enabling us to investigate 
huctuations of the sky brightness at larger angulai_scalns_im to_^M000A_Tlhs_work_js 


comp 


2012 


ementary to recent observational studies fjCoorav et al. 


Zemcov et al. 


2012 b 


Kashlinskv et al. 


201411 in the sense that we also try to hud the existence of the spatial 


huctuations at larger angular scale, but at diherent wavelengths and diherent regions of the 
sky. The huctuation spectra over a wide range of angular scale as well as the spectral shape 
would be important for the interpretation of the nature of the near-infrared background 
radiation. 


This paper is organized as follows. In §2, we present observational data of the NEP 
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deep field, pre-processing of raw image, and other procedures for data reduction. In §3, 
we describe power spectrum analysis to measure fluctuations of the sky brightness, and 
we discuss the possibility of fluctuations by instrumental effects in §4. Correction of 
fluctuation spectrum and error estimation are treated in §5. The contribution of foreground 
components is described in §6. In §7, the result of this study is discussed. We summarize 
our results in the hnal section. 


2. OBSERVATIONS AND DATA REDUCTION 

2.1. Observations 

AKARI is an infrared space telescope developed and launched by Japan Aerospace 


Exploration Agency fJAXA 
described in 


. Deta ils of the mission and it s focal plane instruments are 


Murakami et al. 

O 

O 

) and 

Onaka et al. 

(2007) 


(120071) . AKARI carried out a number 


of wide area surveys in near-, mi d-, and far-infrared wavelength bands. / 


and deep surveys in NEP region flMatsnhara et al. 


2006 


Wada et al. 


2008 


KARI made wic 


Kim et al. 


201211 . 


The survey of the NEP deep held started in 2006 May with all nine wa velength bands of 


Infra-Red Camera (IRC) that is one of the focal plane instruments (see 


Onaka et al. 


2007 


for the detailed description of this instrument) ranging from 2.4 to 24 /rm, and performed 
about 270 pointed observations over one year period, covering nearly a circular area of 0.4 
square degrees (Figured]). Among three near-infrared wavelength bands, we used N2 (2.4 
fim) and N3 (3.2 pm) for this study, as N4 (4.1 pm) had shorter exposure time than the 
other ones. Astronomical Observation Template (AOT) of the NEP deep held is IRC05 
that performed hve sets of exposures (5 images) without dithering for one near-infrared 
wavelength band during one pointed observation. The exposure time for an image is 65.4 
seconds. The held of view (FOV) of an image is about 10' x 10', and the pixe l scale is 1'.'46. 


Further information for observation is available in the IRC data user manual flLorente et al. 
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20081 ). 


2.2. Pre-processing 


We used A KARI IRC iinaging pipeline (version; 070908) for most of the pre-processing 


of each image flLorente et ah 


20081) . However, a few procedures, such as dark subtraction. 


flat held correction, and aspect ratio resampling, were treated separately to improve the 
quality of the image. In the ordinary pipeline, superdark frame that was made by combining 
about hundred dark frames taken during Large Magellanic Cloud (LMC) observations is 
used as a default to subtract dark current. However, ordinary pipeline with superdark 
could not be used to estimate different enhancement of dark current at each pixel in 
the arr ay when the telesc o pe wa s exposed to a large number of cosmic rays. For this 


reason. 


Tsumura fc Wadal (120111) studied a new dark subtraction method that is capable of 


estimating different response of dark current at each pixel even when images were exposed 
to many cosmic rays. 

For the correction of hat held, we made self hat templates using all the NEP deep 
held images taken during the period when stray light was negligible. About 210 images 
in the NEP deep held were used to make self hat template for each band. To make self 
hat template, bright pixels in the image were masked by clipping procedure to remove real 
objects, cosmic rays, and some prominen t artifacts such a s multiplexer bleed (MUXbleed 


hereafter) trails near very bright objects flKim et al.l 120121) . On the other hand, faint pixels 


that were ahected by various kinds of artifacts were treated manually. We examined every 
image and manually masked faint artifacts such as bad pixels, column pulldown, and center 
of saturated objects. Eventually, we obtained self hat templates that show improved feature 
of hat held and describe systematically faint regions in all images. Some stripes, which 
can induce artihcial huctuations toward the horizontal direction of each image, are seen 












after applying the correction for the aspect ratio. To prevent such stripe-like artifacts, we 
adopted the “nearest” option instead of “linear” option for resampling of pixel values. 


2.3. Correction of Instrumental and Time Varying Components 

After pre-processing, additional correction procedures to remove instrumental effect 
and zodiacal light contribution are needed. The hrst additional procedure is related to 
MUXbleed. Observed images of the NEP deep held contain many bright objects, and 
some of them induced MUXbleed. MUXbleed manifests as bright stripes along horizontal 
direction of each image caused by bright objects or strong cosmic rays. MUXbleed stripes 
do not give any signihcant problem in this study because they can be easily removed by 
clipping procedure. However, we found that the pixels in the upper region of MUXbleed 
trail become fainter compared to those in the lower region. In Figure [2l one example of 
an image (pointing ID: 2110163) and average pixel value of background for each row is 
presented to illustrate MUXbleed problem. As shown, there is a large jump in pixel values 
across the MUXbleed trail. Unfortunately, MUXbleed occurs rather frequently in the NEP 
deep imaging data. Among the selected images for this study, about 76% and 78% for 2.4 
and 3.2 /rm, respectively, suffer from the MUXbleed effect. It would be best to exclude 
all pixels whose values are affected by MUXbleed, but we need to use such pixels in order 
not to reject too many pixels. For this reason, we tried to correct for MUXbleed rather 
than abandoning all affected pixels. In Figure El left image and graph show a state before 
correcting for the MUXbleed effect, whereas right parts show a state after the correction. 
Open circle represents average background value of each row, and dotted line in the graph 
shows the position of MUXbleed. We assumed that the effect of the MUXbleed reaches 
the dot-dashed line. Solid line in the left graph is the linear £t of the data points that are 
thought to have been affected, whereas dashed line is the linear fit of data points that are 


regarded not to have been affected by MUXbleed. In other words, dashed line is treated as 
a value of the background of the image when MUXbleed problem does not exist. Correction 
for the MUXbleed problem was done by using the difference between solid line and dashed 
line in the left graph. In Section HI we try to confirm the validity of the result of MUXbleed 
correction by using the area where the effect of MUXbleed is small in the mosaic image. 


The second additional procedure is related to the zodiacal light. In order to make 
an accurate mosaic image, we have to correct for the seasonal variation of the zodiacal 
light. Because the axis of Earth’s orbit is slightly inclined from that of IPD cloud, the 
zodiacal light towards the NEP shows variation over the observational period. To correct 
for such effect, we tried to ht the sky brightness by using sinusoidal function with the 


assumpti on that the spatia' 


latitude flKelsall et al 


1998 


distribution of t he zodiacal light is smooth at high ecliptic 


Pvo et al. 


201211 . The explicit form for the intensity of the 


zodiacal light as a function of the Earth’s heliocentric ecliptic longitude, I, can be written 


I{1) = asin(/ — b) + c, 


( 1 ) 


where parameters a, b, and c represent amplitude, phase angle, and mean sky brightness. 


.especively. wh.ch.,M,b^te^„ed by fitting the obsetved sby bt.gfitness. We used 


MPFITU package flMarkwardt 


20091) for sinusoidal htting, and results are presented in 


Table [H In Figure |3l we show the htting result together with the sky brightness data at 
two near-infrared bands. We subtracted the htted function from each image to correct for 
the time variation of zodiacal light. 


The MUXbleed and time variation of zodiacal light may induce artihcial spatial 
huctuation of sky brightness. Moreover, when we make a mosaic image, sky level of each 
image is adjusted by using overlapped regions between neighboring images. Therefore, two 


%ttp://purl, com/net/mpht 
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additional correction procedures that stated above should be applied before making mosaic 
image. 


2.4. Source Subtraction 


According to 


Takagi et al 


fl2ni2h and 


Kato et al 


(120121) . Point Spread Function (PSF) 


of AKARI IRC images is not circular but elongated along a certain direction due to 
jitters in the satellite pointing, and the direction of elongation varied depending on the 
position in the mosaic image due to variation of position angle of each image. Therefore, 
source subtraction should be done before making a mosaic image. Removing procedure of 
resolved objects in the image consists of two parts: clipping and PSF subtraction. Specific 
procedures for source r emoval are almost the same as those of the previous study for the 


Monitor field fjPaper ll) . but we applied additional procedures for saturated bright objects 


and extended objects to minimize the possible contamination of the relatively faint wing 
parts of such objects. For saturated objects, sky level of an individual object was estimated 
using a mean of pixel values within an annulus whose inner and outer radii are proportional 
to the brightness of the object, and then we additionally applied a circular mask whose 


radius was extended until residual of the object becomes sky level. Fo r extended o 


that were identified by stellarity index in the catalog of optical image flHwang et al 


jjects 


2003 ), 


we also applied additional circular masks with the same ways for saturated objects. 


2.5. Mosaic Images 

AKARI NEP deep survey observations started in 2006 May, but data taken in the 


northern sky 
(IWada et ahl 


durin g May, June, and July 2006 were affected by scattered Earth shine 
20071 ). For this reason, we excluded the data taken between April and August 
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to avoid the Earth shine problem. We selected 86 pointed observations for two bands taken 
from 2006 September to 2007 March. We further rejected additional images to obtain a 
reliable result of this study. First, images containing the planetary nebula NGC 6543 (the 
bright region in Figure [T]) were rejected because it is not easy to subtract the contribution 
of a bright and irregular object. Second, we did not use images in the inner annulus but 
used only those in the outer annulus because of the difficulties in making a reliable mosaic 
image that includes both inner and outer annuli. To obtain a stable and reliable mosaic 
image, adjacent images should share sufficient FOV. As shown in Figure [U however, the 
overlapping region between inner and outer annuli is too small for that. Third, bad images 
that show poor stability during observation or various artifacts including cosmic rays were 
also rejected. Eventually we used 75 and 68 images for 2.4 and 3.2 /im, respectively, out 
of 34 pointed observations. These images after source subtraction with pre-processing and 
two additional corrections were combined to make mosaic images for two bands by using 
the software Montage^. Montage adjusts sky level of each image by using overlapping 
regions among adjacent images during mosaicking. Because sky level of an image is usually 
not exactly consistent with that of the neighboring image in the overlapping region, sky 
adjustment procedure while making mosaic images further reduces artihcial fluctuations. 

In order to guarantee a stable mosaicking procedure, sufficient overlapping regions 
among adjacent images are necessary. For this reason, although source subtraction 
procedure was performed before mosaicking, further masking for saturated bright objects 
and extended objects was conducted in the mosaic image to leave sufficient overlapping 
regions for mosaicking. Moreover, instrumental and observational artihcial components, 
such as bad pixels, column pulldown, center of the saturated objects, and MUXbleed, were 
removed in advance before mosaicking because they hinder sky matching process because 


^http://montage.ipac.caltech.edu/ for more information about this software. 
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of their different pixel values and shapes in the adjacent images. After mosaicking, pixels 
that remain at least in two pointed observations are selected for the study (the top panel of 
Figure HD . Because pixels that are survived only in a single pointed observation do not have 
a chance to suppress the contributions of artihcial components, they are more likely to be 
affected by artifacts such as leakage of stray light (see Section |4] for details). In addition, 
it is hard to recognize whether the boundary region of an image is affected by objects just 
outside of the image. In order to avoid such contamination, we selected pixels that are also 
survived at least in an additional single adjacent pointed observation. 

The small area denoted by a dashed line box in Figure 0] represents the region where 
the effect for MUXbleed is small compared to the other regions because this area does not 
have any bright object. We use this small area to conhrm the validity of the correction for 
MUXbleed problem. We can measure fluctuations of the sky brightness by using mosaic 
images. Basic parameters of the mosaic image of the NEP deep held after the correction 
procedures are listed in Tabled In the bottom panel of Figure IH smoothed mosaic images 
made by averaging of neighbor pixels within a circular shape hlter with a diameter of 100" 
are presented to show spatial structure more clearly. Although, spatial structure is seen 
in the smoothed mosaic images, its pattern appears to be somewhat different for each 
band. This could imply a possibility of the existence of artihcial structure. More careful 
assessment of such phenomenon is discussed later in Section 01 


3. POWER SPECTRUM ANALYSIS 


We used power spectrum analysis to measure 


held following the previous works (refer to 


spatial huctuation of the NEP deep 


Paper 1 


for the more detailed description of 


power spectrum analysis). The top panel of Figure [5] represents the huctuation spectra, 
[g^F2(9)/27r]^/^, for the mosaic images (hlled circles) and subset analysis (open circles). 
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Filled circles contain contributions of instrumental noise and unresolved faint objects as well 
as background component. We estimated the contribution of each component separately 
to extract background component. Subset analysis is intended to estimate fluctuation 
spectrum of instrumental noise in the mosaic image and to conhrm whether observed 
structure in the mosaic image is of celestial or cosmic origin. To perform subset analysis, the 
whole images were divided into two subsets composed of alternating pointed observations in 
sequence from the left side of the mosaic image, and a difference of two subsets was taken. 
The difference of two subsets will erase persistent features in the sky, and leave instrumental 
noise that consists of random noise and artihcial feature in the mosaic image. As shown in 
the top panel of Figure [5l the subset analysis reveals signihcantly smaller fluctuation than 
that of the mosaic image. However, the result of subset analysis shows larger fluctuation 
than that of random noise at large angular scale. As discussed later in Section 01 this 
implies that there could remain artihcial structure in the mosaic image. A difference of 
huctuation spectra between the subset analysis and random noise represents the amount of 
the contamination by the possible artihcial structure. At angular scales larger than 100 ", 
a mean power spectrum, ^ 2(^)5 due to artihcial structure reaches up to about 10 % and 
12% of that of the mosaic image for 2.4 and 3.2 /rm, respectively. Fluctuation spectra of 
the mosaic images after quadratically subtracting the results of the subset analysis are 
presented in the bottom panel of Figure [5l 


4. CONTRIBUTION OF INSTRUMENTAL EFFECTS 

One of the possible origins of the spatial huctuations of the infrared background is the 
instrumental ehect. MUXbleed ehect could induce artihcial huctuation. In order to check 
any possible artihcial huctuations due to insufficient correction for MUXbleed, we used the 
small area (denoted by dashed line box in Figure 0]) in the mosaic image as a representative 
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region where the effect of MUXbleed problem is negligible even before MUXbleed correction. 
In Figure El two new fluctuation spectra (solid lines and asterisks) are shown along with 
the spectra derived in the previous section (open circles). Solid lines represent fluctuation 
spectra of the small areas where the effect of MUXbleed is negligible compared to other 
regions in the mosaic image. In addition, we made another kind of mosaic image for each 
band composed of images before the correction of MUXbleed. Their fluctuation spectra are 
presented as asterisks. Because the effect of MUXbleed is not confined in a single image but 
propagates to all regions during the mosaicking procedure, asterisks show larger fluctuations 
than open circles at angular scales larger than 100". Moreover, open circles have similar 
fluctuation power to that of solid lines at angular scales up to a few hundreds arcseconds, 
implying that significant amount of MUXbleed problem has been removed by the correction 
procedure shown in Figure El We can see that the whole mosaic image after MUXbleed 
correction shows similar fluctuation power to the small area where the effect of MUXbleed 
is negligible. We conclude that the correction for MUXbleed problem was effective. 

Meanwhile, as shown in the top panel of Figure [5l the result of subset analysis gives 
larger fluctuation than that of random noise at large angular scales. This means that 
subsequent pointed observations toward the same field show different structures with each 
other, which may have been originated from temporary instrumental effects. Although 
contributions of artificial components are estimated and subtracted from the final result by 
subset analysis, we attempted to identify the cause of large fluctuation in subset analysis 
to ensure the result of this study. We made two subsets that contain alternating images 
(i.e., alternating exposures in a pointed observation), rather than alternating pointed 
observations, in sequence from the left side of the mosaic image. As presented with asterisks 
in Figure [71 the result of new subset analysis shows smaller fluctuation than the previous 
one and becomes closer to the fluctuation of random noise shown as a dotted line at large 
angles. This implies that there exist artificial structures in each pointed observation, which 
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persist at least during a period of one pointed observation. In the new subset analysis, an 
image in one subset shares the same FOV with its counterpart in the other subset due to 
no dithering between them. Therefore although there are artihcial structures in each image, 
they are canceled out each other by taking a difference of two subsets if artihcial structures 
are unchanged during a period of one pointed observation. As a result, new subset analysis 
shows nearly random-noise-like huctuation spectrum. 


There are several candidates that can induce artihcial structures in each pointed 
observation. First, errors occurred during pre-processing of raw imaging data can be 
a source of the artihcial structure. Howe ver, beca use we used the same pre-processing 


method with that of the previous study (jPaper ll ) and the result of subset analysis in 


the previous study shows random-noise-like behavior, the problems in pre-processing are 
not likely to be the cause. Another possible origin is the error that arises during source 
subtraction procedures by using PSF. Because PSF of AKARI IRC image is not symmetric 
but elongated and the direction of PSF elongation in the mosaic image varies depending 
on the position angle of the pointed observation, there is a possibility of distinct residuals 
in each pointed observation after PSF subtraction. If so, the ehect of PSF subtraction 
may produce artihcial huctuation in subset analysis at small angular scales of a few tens 
arcseconds as well as at large angular scales. As shown in Figure [71 the large huctuation 
in subset analysis mainly appears at angular scales larger than 100", whereas the result of 
subset analysis follows the behavior of random noise at small angular scales. Therefore, the 
error related with PSF subtraction procedure is not likely to be the main reason for the 
large huctuation at large angular scale in subset analysis. 


The remaining probable cause is temporary instrumental ehects such as stray light that 
remain as distinct structures in each pointed observation. We paid particular attention to 
the fact that stray light can change during a period of subsequent pointed observations. To 
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test the stray light as a possible origin of the artificial fluctuation, we made more detailed 
analysis of the imaging data of the Monitor field. We found 13 pointed observations of the 
Monitor field that were carried out during the same period, and 13 images were made by 
stacking 3 images of N2 band in each pointed observation separately. After that, 13 images 
were subtracted from the entire stacked image of N2 band in the Monitor field separately. 
We present result in Figure [HI which shows distinct structure in each pointed observation. 
They are artificial structures because real structure is removed by subtraction of the entire 
stacked image of the Monitor field. The last image in Figure [8] was made by stacking 13 
images, which shows negligible structure and indicates that artificial structures in each 
pointed observation have random nature compared to those in other pointed observations. 
The test with the Monitor held data has some advantages in identifying possible candidates 
that could lead to the large huctuation in subset analysis. For example, because they do not 
contain MUXbleed problem a nd subset analysis of the Monitor held shows similar behavior 


with random noise fjPaper ll ). we can exclude MUXbleed problem and errors occurred 


during pre-processing from the possible candidates. In addition, because only clipping 
procedure was used to remove objects in this test, we can also exclude errors occurred 
during PSF subtraction procedure as a possible candidates. 


We made same test for N3 band to examine the behavior of artihcial structures in the 
Monitor held. Note that AOT of the Monitor held is IRC03 that is diherent with that of the 
NEP deep held. IRC03 performed three sets of images for three near-infrared wavelength 
bands of the IRC alternately during one pointed observation, and hlter wheel change or 
dithering was applied between exposures. We also found artihcial structures in N3 band and 
obtained linear Pearson correlation coefhcients of 0.112(mean)±0.158(standard deviation) 
with those of N2 band in each pointed observation. Evidently there is no statistically 
signihcant correlation between them even in each pointed observation. That is to say, 
unlike artihcial structures in the NEP deep held, those in the Monitor held vary frequently 
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even during a period of one pointed observation. A probable candidate that can explain 
both the existence of artificial structures in two fields and the difference of their behaviors 
is a contribution of stray light that depends on the position of filter wheel. Although we 
used imaging data that were observed during a period of small stray light, we could not 
exclude some leakage of stray light that has contaminated the images. In the case of the 
NEP deep held, because there is no operation of hlter wheel during a pointed observation, 
the pattern of the stray light remained nearly the same for different images taken during 
one pointed observation. On the other hand, in the case of the Monitor held, because 
stray light could have been scattered diherently due to operation of hlter wheel during a 
pointed observation, we can expect diherent shape of artihcial structure in each image. 
Therefore, distinct artihcial pattern may remain in each pointed observation of the NEP 
deep held more clearly than that of the Monitor held. Moreover, such artifacts are rarely 
suppressed during mosaicking due to small number of stacked images per held. To evaluate 
the contribution of artihcial structures quantitatively, we estimated huctuation spectra of 
the images in Figure [HI The results for 13 pointed observations are presented as solid lines 
in Figure ini and a mean of them is presented as hlled circles, which show larger huctuation 
than dotted line at large angular scales. Fluctuation spectrum of the last image denoted by 
open circles shows small deviation from random noise and smaller overall huctuation power 
than those of other images. 

The amount of contribution of artihcial structure in the mosaic image of the NEP deep 
held can be estimated by quadratically subtracting the huctuation spectrum of random 
noise from that of subset analysis in the top panel of Figure [5], which becomes 0.098 and 
0.053 nW m~^ sr“^ for 2.4 and 3.2 /xm, respectively, as a mean at angles between 100" 
and 350". We can expect such quantitative contribution of artihcial structure of the NEP 
deep held by converting that of the Monitor held presented in Figure [9l A diherence 
between hlled circle and dotted line represents huctuation due to artihcial structure in 
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each pointed observation of the Monitor held. For the case of N2 band, a mean difference 
at angular scales larger than 100" is about 0.092 nW sr“^. In a pointed observation 
of the Monitor held, an amplitude of artihcial structure becomes \/3 times weaker as 3 
images with diherent shape of artihcial structure are stacked with dithering. However, such 
a reduction would not occur for the NEP deep held because images with similar shape of 
artihcial structure are stacked without dithering. As a result, in the NEP deep held, we can 
expect that the contribution of artihcial structure in each pointed observation is \/3 times 
larger than that of the Monitor held. Finally, because the mosaic image for N2 band in the 
NEP deep held is covered by 3.4 pointed observations per held on average, the contribution 
of artihcial structure becomes ^/3A times weaker in the mosaic image, resulting in 0.086 
nW m“^ sr“^. Similarly we estimated the huctuation due to artihcial structure to be 0.049 
nW m“^ sr“^ for N3 band in the NEP deep held. Such result corresponds to about 90% of 
the residual huctuation of subset analysis over random noise shown in the upper panel of 
Figure |5l Therefore we conclude that the contribution of temporary instrumental ehects, 
most likely stray light, may be the most prominent cause for the large huctuation in subset 
analysis. Although this test is valid only within angular scales of a few hundreds arcseconds, 
it can give a hint for the large huctuation of subset analysis at angular scales larger than a 
few hundreds arcseconds. 

Although we examined artihcial components in the mosaic image and their possible 
origin, it is difficult to correct them more accurately due to their random nature in each 
pointed observation. To reduce further the contribution of artifacts with random nature, 
we need larger number of stacked images per held. 
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5. CORRECTION AND ERROR ESTIMATION 
5.1. Corrections of the Fluctuation Spectrum 


Fluctuation spectra in the bottom panel of Figure [5] were further corrected to remove 
instrumental and systematic effects. Corrections of the fluctuation spectra were conducted 
as follows. The power spectrum after the corrections can be expressed as 


P2{q) 


M-\q)P,{q) 

T\q)B{q) 


( 2 ) 


where M{q) is a mode coupling matrix, T{q) is a map-making transfer function, B{q) is a 
beam transfer function, and ^ 2 ( 9 ) represents power spectrum before the corrections. 


Observed fluctuation is affected by the presence of mask in the mosaic image. In 


order to correct the observed fluctu ation, we examined an effect of the mas. 


i on each 


data point separately by simulation flCoorav et al. 


2012b; 


Zemcov et al 


20141) . We made 


a simulated image that has a unit value of fluctuation power at a certain angular scale 
but no fluctuations at other angular scales. By adopting mask of the mosaic image on 
the simulated image, we can estimate how the mask converts the unit fluctuation power 
at a certain angular scale into those of other angular scales. Simulations were iterated 
100 times for each data point by adjusting phase of a complex number in Fourier space, 
and a result was expressed as a matrix, M{q), whose each row corresponds to the result 
of simulation for each data point. The mode coupling matrix can be used to estimate 
the effect of mask on the fluctuation, whereas its inverse, M~^{q), can be used to recover 
the intrinsic fluctuation. To check usefulness of the mode coupling matrix, we introduced 
a seed fluctuation that has similar fluctuation power to the observation and produced 
500 simulated images with identical fluctuation spectrum to the seed fluctuation. Each 
simulated image corresponds to a random realization of the seed fluctuation by adjusting 
phase in Fourier space. Fluctuation spectra of the simulated images with mask of the mosaic 
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image were recovered by the inverse of the mode coupling matrix. As shown in Figure fTOl 
mean of recovered fluctuations (solid line) is well matched with the seed fluctuation (open 
square). By using simulated mask pattern, we conhrmed that an error of the recovered 
fluctuation is related with the portion of the mask. Because FOV of the NEP deep held is 
not rectangular, outside of the FOV is also treated as a mask in the simulation. Therefore 
such large portion of the mask induces large error of the recovered huctuation. Figure dO] 
shows that the mode coupling matrix is efficient in estimating the mean effect of the mask 
for multiple realizations. Because there is no information about an intrinsic huctuation 
spectrum of the NEP deep held, we do not know how the intrinsic huctuation spectrum of 
the NEP deep held is to be converted by mask of the mosaic image. Therefore, we adopted 
the mean ehect of the mask for the correction. Two components are considered for the 
error of the mode coupling matrix. One is the systematic error that can be measured by 
diherence between seed huctuation and recovered huctuation in Fig. 10. The other is the 
statistical error as measured by the standard deviation of the recovered huctuations for 
many diherent realizations. Two components are added quadratically, and it becomes the 
most dominant component of the error budget at angular scales larger than 300". However, 
at that angular scale, mean errors of the mode coupling matrix for 2.4 and 3.2 ^nm, 0.086 
and 0.054 nW m“^ sr“^, are smaller than observed huctuations for 2.4 and 3.2 ^m, 0.177 
and 0.129 nW m“^ sr“^. 

Fluctuation power at small angular scale can be reduced due to the ehect of PSF. 
Although we used analytic form for the PSF to subtract wing parts of detected objects, 
it is difficult to dehne asymmetric central parts of objects by using the analytic function. 
Therefore, to correct for the ehect of PSF on the huctuation, we derived PSF of the NEP 
deep held by stacking point objects. We stacked 127 and 85 point objects for 2.4 and 3.2 
/rm, respectively, to obtain PSFs and their power spectra of these bands. The beam transfer 
function in Figure ITT] represents power spectrum of PSF that is normalized to be 1 at large 
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angular scale, and error bar represents a variation of PSF in pointed observations. Because 
PSF of one pointed observation is a little different from those of other pointed observations 
due to variation of jitters in the satellite pointing, we modeled PSF and estimated its power 
spectrum separately for each pointed observation. Error bar in Figure fTTl represents their 
standard deviation. To correct for the contribution of PSF, observed power spectrum was 
divided by the beam transfer functio n. The error was included in t he hnal error budget after 


2 


1992 


, page 43). At angular 


applying error propagation formula fiBevington fc Robinsor 
scales larger than 100", the error related to the beam transfer function becomes smaller 
than 0.002 and 0.001 nW m“^ sr“^ for 2.4 and 3.2 pm, respectively, which are smaller than 
observed fluctuations by two orders of magnitude and negligible compared to other error 
components. 


In order to test a rel iability of the 


transfer function, T{q) fjCoorav et al 


mosaic king procedure, one usually uses a map-making 


2012b|) . First, we generated an image of the sky 


that is large enough to cover the analyzed portion of the NEP deep held in this study (36' 
X 25' and 35' x 23' for 2.4 and 3.2 pm, respectively) with huctuation power similar to 
the observed one in order to mimic the real sky. Individual frames identical to the size of 
AKARI IRC frame have been extracted from the simulated image. We then assigned the 
same noise and sky levels, and applied the same masking patterns with observed frames 
to all extracted frames. An intensity of the random noise in the simulated image was 
determined by the result of subset analysis at small angular scale, and sky level of each 
simulated image was determined using the value in Figure El Finally, they were combined 
through the same mosaicking procedure using Montage as we did for the observed images. 
Simulations were repeated for 1000 times, and original huctuation spectra were compared 
with those of remade mosaic images. The map-making transfer function in Figure [12] 
represents mean of the ratio of remade huctuation spectra to original ones, and error bar 
represents a standard deviation of 1000 simulations. First of all, at large angular scale. 
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we were able to reproduce the original fluctuation for 2.4 /rm, while we found fluctuation 
larger by about 10% at the largest angular scale for 3.2 pm. One possible cause of such 
a difference is difference in overlapping regions between 2.4 and 3.2 pm: there are more 
overlapping region in 2.4 pm than in 3.2 pm. However, this is mainly due to the order of 
data reduction procedure in this study. Because we masked some pixels before making 
the mosaic image to remove objects, overlapping region between adjacent images becomes 
small. It may be a main reason that for additional artificial fluctuations at large angular 
scale for 3.2 pm. To confirm this, we obtained the map-making transfer function for 3.2 pm 
without masking before making the mosaic image. At large angular scale, we found that 
the additional fluctuation is reduced to about 2% of the original fluctuation for 3.2 pm. 
Also the standard deviation of 1000 simulations decreased by about one third compared 
to the result in Figure [121 On the other hand, significant suppression in the fluctuation 
spectrum is apparent at small angular scale due to smoothing effect by resampling of pixel 
values during mosaicking procedure. We found that such suppression is mainly induced by 
instrumental random noise rather than celestial objects at angular scales smaller than 10". 
Therefore, we did not use the fluctuation spectrum at angular scales smaller than 10" for 
further analysis since the measured fluctuations at such small scales are dominated by the 
shot noise and he beyond our interests. To correct for the effect of the map-making transfer 
function, observed fluctuation spectrum was divided by the map-making transfer function. 
The error was included in the final error budget after applying the error propagation 
formula. At angular scales larger than 100", mean of the error related to the map-making 
transfer function is 0.009 and 0.006 nW m“^ sr“^ for 2.4 and 3.2 pm, respectively, which 
are smaller than observed fluctuations by an order of magnitude. 
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5.2. Systematic Errors 

In the previous section, errors related with the mode coupling matrix, the beam 
transfer function, and the map-making transfer function were described. In addition, flat 
held error and an error due to residual of detected objects were investigated as possible 
systematic errors because their contributions to the observed huctuation spectrum may be 
signihcant. Even a small amount of hat held error in an image can be accumulated in the 
mosaic image during mosaicking procedure. Therefore, artihcial huctuation due to hat held 
error can ahect the hnal huctuation spectrum. In order to obtain hat held error, we stacked 
190 and 220 images for 2.4 and 3.2 ;um, respectively, with (x, y) coordinate after masking 
detected objects. After that, the obtained hat pattern is assigned to individual frames for 
each band and mask pattern and sky level were also assigned to them. Finally, they were 
combined through the same mosaicking procedure using Montage as we did for the observed 
images, and huctuation spectra were estimated. The result of this simulation can be treated 
as an upper limit due to hat held error because simulated image also contains celestial 
huctuation and huctuations due to other artifacts. As shown in Figure [131 huctuation 
spectra of the simulated images are smaller than observed huctuations at all angular scales, 
indicating that hat held error is not signihcant and not severely accumulated in the mosaic 
image during mosaicking procedure. 

Residuals after masking and subtraction of detected objects can also give some 
contributions to the hnal result. If there remain residuals of detected objects, their 
contributions are most signihcant just outside of the masked region of the objects. For this 
reason, we examined residual within 5 pixel layers from the masked region by comparing 
sky level of each detected object. After that, to make a simulated image with a size of the 
mosaic image, mean of residual of each detected object is assigned to that region, while 
zero value is assigned to other regions. Instead of real residual, mean of residual is assigned 
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to suppress the contribution of noise and undetected faint object in that region. As shown 
in Figure [T31 the contribution from the residual of detected objects is also smaller than 
observation especially at large angular scale, implying that the observed excess fluctuation 
at large angular scale is not affected by the residual of the detected objects. 


5.3. Corrected Fluctuation Spectrum 


Fluctuation spectra after corrections by the mode coupling matrix, the beam transfer 
function, and the map-making transfer function are presented in Figure fT3l As shown, the 
mosaic images have excess fluctuations, i.e., residual fluctuations over the shot noise, at 
large angular scales. Errors related with corrections and two systematic errors are also 
plotted separately along with la error of Fourier elements of each bin. Error bars of the 
hlled circles are quadratic sum of them. The effect of shot noise originating from unresolved 
faint galaxies was estimated by using simulated images of the same size with the mosaic 
image containing objects fai nter than the limitin g magnitude. We used the power-law 
distribution, N{m) oc (IMaihara et al.ll200ll) . of the number of artihcial objects of 

magnitude, m, in unit magnitude interval for the generation of simulated images. The 
normalization was done based on the number of detected objects in the mosaic imag e. 


ARTDATA package of the Image Reduction and Analysis Facility (IRAF; 


Tody 


19861) was 


used to generate a list of randomly distributed objects and to assign them on the simulation 
image. Before the estimation of fluctuation spectrum, the mask pattern of the mosaic 
image was applied to the simulation image. We made 100 simulated images for each band 
and estimated mean of 100 fluctuation spectra. The overall amplitude of fluctuation is 
determined by the limiting magnitude. The observed fluctuation spectra at small angular 
scale well match the shot noise with the limiting magnitude (AB) of 21.4 and 21.8 for 2.4 
and 3.2 nm, respectively. 
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6. CONTRIBUTION OF FOREGROUND SOURCES 
6.1. Zodiacal Light 


Because zodiacal light is the most signihcant contributor to brightness of the near- 
infrared sky, it is a probable candidate to explain the excess fluctuation of the mosaic image. 


Pvo et al. 


(120 12|) investigated spatial fluctuation of mid-infrared sky by using the Monitor 


held data set of AKARI and found very small residual huctuation over photon noise. Since 
we expect the huctuation due to zodiacal light does not depend on the wavelength, we 
adopted their results of ~0.03% of the sky brightness at 18 ^m band as an upper limit of 
huctuations due to zodiacal light in the near-infrared sky at angular scales larger than 100". 
If we estimate huctuations by using the results in the mid-infrared sky, the upper limits of 
huctuations due to zodiacal light become 0.034 and 0.021 nW m“^ sr“^ at 2.4 and 3.2 pm, 
respectively. Due to the motion of the Earth, IPD column along the line of sight toward 
the NEP continuously changes during observation period. The huctuation due to zodiacal 
light could have been additionally reduced during mosaicking because subsequent pointed 
observations toward the same held might detect diherent shape of spatial huctuations in 
the zodiacal light. The estimated upper limits due to zodiacal light are clearly smaller than 
the observed excess huctuations. 


6.2. Galactic and Extragalactic Components 


Dihuse Galactic light (DGL), starlight scattered by the interstellar dust, can also be a 
candidate of the excess huctuation. Because DGL traces the distribution of the interstellar 
dust, we examined far-infrared 90 pm band that was taken by Far-Infrared SurveyoJ^ (FIS; 


^FIS images are available at http://darts.jaxa.jp/ir/akari 
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Kawada et al 


2007D on board AKARI. Slow scan observation with a scan speed of 15 


arcseconds per second and a reset interval of 2 seconds were carried ont to cover the NEP 
eep region. We pro cessed the raw pointed data with the official FIS slow-scan toolkit 


(IVerdugo et ah 


20071) . The data rednction tool removed bad data and cosmic rays and 


corrected for the variation of the detector responsivity. After the dark current subtraction 
and flat holding, undetected glitches and other artifacts were removed. The uncer tainties 


of hux calibration is no more than 20% for the 90 fim band flKawada et al 


20071). The 


processed images were combined to make a mosaic image using the software Montage. 


To remove the contribution of discrete objects in the far-infrared mosaic image, we 
masked pixels whose values exceed three times of the standard deviation from mean pixel 
value, and one layer of pixel around the masked region was additionally masked. About 
5% of pixels were masked in this procedure. To estimate the correlation between images 
with different pixel resolutions, the angular resolution of the near-infrared mosaic image 
was adjusted to that of the far-infrared mosaic image, and then a mask of the near-infrared 
mosaic image was applied on the far-infrared image. The left panel of Figure 1141 shows that 
there is no clear correlation between near- and far-infrared images. However, as discussed 
earlier, because there are some artihcial structures in the near-infrared mosaic images, this 
method has some limited implications. For this reason, we made another examination. The 
right panel of Figure [H] shows fluctuation spectra of the far-infrared image at both the 
Monitor held and the NEP deep held. We focused on the huctuations at angular scales 
larger than 100" where the shot noise becomes small. The far-infrared images of both 
helds have similar huctuation powers at angular scales up to 300", and the NEP deep held 
does not show any signihcant increase of huctuation even at larger angles. Such a result 
represents that the NEP deep held has similar contribution of DGL to the Monitor held. 
Therefore, we can say that the contribution of DGL is negligible at the NEP deep held by 
using the result in the right panel of Figure [TH 
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The clustering of unresolved galaxies fainter than the limiting magnitude of the mosaic 
image can also be a candidate for the excess fluctuations. Like the previous study with 


the Monitor held, the feature of obser ved excess huctua' 
shows inconsistency with the study by 


Chary et ah 


ion with blue spectrum also 


(120081) who proposed that a fraction of 


huctuations in the near-infrared background could be induced by red dwarf galaxies. We 
used the deep imaging observations at Kg band (Ac: 2.146 fim) with Canada-France-Hawaii 
elescope (CFH T) Wide-held InfraRed Camera (WIRCamJ^ for the NEP deep held 


flOi et al. 


20141) for the estimation of the contribution of the clustering of unresolved faint 


galaxies to the excess hu ctuations. To extract sources in Kg band image, we used SExtractor 


(IBertin fc Arnouts 


1996 ) with the following parameters: DETECT THRESH=2, DETECT 


MINAREA=3, BACKSIZE=64, BACK FILTERSIZE=3. Detection limit was estimated 
by using simulated sources and the same parameter sets of SExtractor to recover them. 
The 50% detection limit is about 23.5 AB mag that is deeper than the limiting magnitude 
of N2 band and is sufficiently sensitive to detect main contribution from the clustering of 
unresolved faint galaxies. The simulation was performed using the detected objects at Kg 
band in the area (dotted line box in Figured]) that covers most regions used in this study. 
We prepared a simulation image having the same size with the area with dotted line box 
and assigned detected objects at the same position on the simulation image. When we make 
a simulated image with objects fainter than 21.5 AB mag at Kg band, huctuation power of 
the simulation at small angular scale becomes the same as that of N2 band mosaic image or 
shot noise. A result of the simulation is presented as a solid line in Figure dS] which shows 
larger fluctuation than shot noise at large angular scales. We obtained 0.059 nW m~^ sr“^ 
as a mean contribution of the clustering of unresolved faint galaxies at angular scales larger 
than 300", which is 2.8 times smaller than the observed excess fluctuation at the same 


“^http: //www.cfht.hawaii.edu/Instruments/Imaging/WIRCam 




















27 


angular scale. Hence we conclude that the contribution of the unresolved faint low-redshift 
galaxies is not enough to explain the observed excess fluctuations. 


Sullivan et al. 


(120071) also investigated angular power spec tra of res olved sources at J, 


Ks, and L bands. The previous study with the Monitor held fjPaper ll ) used their model 
to estimate the contribution of unresolved faint low-redshift galaxies and showed that the 
estimated contr i butio n is small compared to the observed huctuation at 2.4 /rm. Recently, 


Helgason et al. 


(120121) examined huctuations in the infrared background from known galaxy 


populations by using large sets of luminosity functions. They applied their model to 
actual measurements and concluded that the known galaxy population is not likely to be 
a candidate to explain the observed huctuations at arcminute scales. Although it is hard 
to compare their results with ours quantitatively because our simulation only used the 
detected objects at Ks band without completeness correction at the faint-end of luminosity 
function, every study predicts that the contribution of unresolved faint low-redshift galaxies 
to the huctuation power is small compared to the observational measurements. 


DISCUSSION 


By using observational data, we found excess huctuation in the near-infrared 
background at angular scales up to 1000". We estimated a mean of the excess huctuations of 
0.26 and 0.11 nW m“^ sr“^ for 2.4 and 3.2 pm bands, respectively, at angular scales larger 
than 100". In Figure [161 we plotted the results (open squares) and those of the Monitor held 
(open circles). In addition to the similar level of excess huctuations, the color of huctuation 
components of the NEP deep held is found to be similar to that of the Monitor he ld: both 


of them show blue spectra, as hrst observed by the IRTS (IMatsumoto et al. 
excess sky brightness. 


2005 ) in the 
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By using the result of the NEP deep held, we have extended huctuation analysis of the 
Monitor held at angular scales up to 1000". In Figure [T71 the result of the Monitor held 
(hlled circles) is combined with that of the NEP deep held (open squares). The huctuation 
spectrum of the NEP deep held connects smoothly with that of the Monitor held and 
extends without showing any signihcant variation of huctuation power at angular scales 
larger than a few hundreds arcseconds. Fo r comparison with the result of AKARI, that of 


Spitzer at 3.6 /rm flKashlinskv et al. 


2 OI 2 I) is plotted, which is adjusted by the assumption 


that the spectrum of huctuating components at the near-infrared wavelength is close to 
Rayleigh-Jeans spectrum (~A“^). It is seen that the average excess huctuation of AKARI 
is consistent with the result of Spitzer at angular scales between 100" and 1000". 

Although the origin of such excess huctuation is not yet well understood, we could 
exclude zodiacal light, DGL, and unresolved faint galaxies at low-redshift. One of the 
remaining possibilities is the contribution from the very early objects at z > 10. Theoretical 
studies have been performed to estimate huctuation due to the clustering of high-redshift 
objects. Their results show that huctuations of high-redshift objects are diherent from 
those of low-redshift objects. Our result can be explained from this point of view. Some 


theoretical studies wit 


1 biased halo mode 


feature at around 10' (jCoorav et al. 


2004; 


show the huctuation powers with a turnover 


Kashlinskv et al. 


predict gradually decreasin 


; huctuation powers t o ward 


nonlinear biasing of halos flFernandez et al 


2010 


2012 


2OO4J) . while some other studies 
arge angular scales by considering 


Coorav et al. 


2012al) . The result 


of this study prefers the latter case, as shown in Figure [IZl because our result does 
not show any signihcant turnover feature at up to around 10'. However, it should be 
interpreted carefully because there are only a small number of Fourier elements in the 
last few data points. Meanwhile, although TeV 7 -ray observation prefers lower level of 


background brightness than observational measureme nts in the near-infrare c 


explain absorption feature in 7 -ray spectra of blazars flAharonian et al. 


2006; 


bands to 


Gilmore et al. 
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20121: iMever et al.l 120121: IH.E.S.S. Collaboration! 120131) , it is not a stringent constraint 
because alternative explanation of 7 -ray spectra of blazars is possible using sec ondary TeV 
7 -ray photons produced by in teraction between cosmic rays of blazars and EBL flEssev et al. 


2010 


Essev fc Kusenko 


20101 ) . 


However, theoretical studies p redicted smaller flu ctuations than observations by about 


an order of magnitude. Recently, 


Yue et al. 


fl2013af) investigated the contribution of 


high-redshift galaxies at z > 5 to the near-infrared background by using simulations with 
an analytical model. They predicted that fluctuation amplitude due to 7 > 5 galaxies 
would be less than 0.01 nW m“^ sr“^ at a few hundreds arcseconds in the wavelength range 
1.0~4.5 /rm, which is 10 times smaller than the ex cess fluctuations in this study and rece nt 


obser vational studies using Spitzer imaging data flCoorav et al. 


2012 b 


Kashlinskv et al. 


20121 ). Therefore, further studies for both theory and observation are required to identify 


the origin of excess fluctuation in the near-infrared background. 


Coorav et al. 


(j2012bl) suggested that the extended emission due to IHL can be a 


candidate to explain observed fluctuation in the near-infrared background. Although it is an 
interesting idea, it may be challenging to verify observationally. As discussed in their work, 
simple masking is not an efficient method to remove diffuse components of faint galaxies. 
To identify extended components around faint galaxies and to subtract them accurately, we 
need imaging data with high sensitivity and a sophisticated method for modeling them. 

Recently, the rocket-borne instrument, GIBER, measur ed the fluctuation o f the sky 


brightness up to 1° angular scale at 1.1 and 1.6 ^m bands fjZemcov et al. 


20141) . Since 


there exist fluctuations at shorter wavelength than 2 pm, IHL model is preferred than the 
contribution from Population HI stars in the early universe. They also found that the result 
of GIBER shows good correlation with the result of Spitzer and fluctuation at 1.1 pm shows 
some departure from the Rayleigh-Jeans spectrum. 
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Yue et al.l fj2013b|) suggested a contribution of direct collapse black holes (DCBHs) 
at 2 ; > 12 as a possible candidate for near-infrared fluctuations. They claimed that the 
DCBH model can explain observed fluctuations at 3.6 and 4.5 /mi and observed cosmic 
X-ray background/near-infrared background cross-correlation. However, DCBH model 
i s inconsistent with substantial fluctuations even around 1 /im observed by the GIBER 


fjZemcov et ah 


201411 . 


To investigate various aspects of the excess fluctuation, more observational studies are 
needed by using individual imaging data with various FOV, angular scales, and wavelength 
bands. As one of the scheduled programs for the st udy of the near-infrare d background. 


the Multipurpose Infrared Imaging System (MIRIS; 


Han et ah 


2010 


20141) . launched in 


Nov. 2013, is observing a large area of 10° x 10° centered at NEP for the study of the sky 
brightness and its fluctuations at 1.1 and 1.6 /im bands. It will provide new information 
about the sky brightness in the near-infrared bands at angles up to a few degrees and 
improve our understanding of the origin of the cosmic infrared background. 


8. SUMMARY 

We studied spatial fluctuations of the sky brightness using the NEP deep survey data 
obtained with AKARI IRC at two wavelength bands, 2.4 and 3.2 //m. Compared with the 
Monitor held of AKARI used in our previous study, the NEP deep held has larger areal 
coverage, enabling us to extend huctuation analysis of the sky brightness at larger angular 
scales up to 1000". We produced mosaic images for two bands and estimated spatial 
huctuations using power spectrum analysis. We found that there is a residual huctuation 
over the estimated shot noise even at larger angles than the angular scale of the Monitor 
held, which can hardly be explained by any known foreground components. 
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By using excess fluctuation of the NEP deep field, we confirm a blue stellar-like 
spectrum of excess fluctuation that was found in the previous study with the Monitor held 
and hud the blue spectrum is still recognized even at larger angular scales. The excess 
huctuation spectra of the NEP deep held extends smoothly that of the Monitor held at 
angular scales of a few hundreds arcseconds, and combination of the two results does not 
show any signihcant variation of huctuation power at angular scales between a few hundreds 
and ~1000". We show that the result of this study is consistent with the study of Spitzer 
at 3.6 pLvu. Such results could provide useful constraints for the study of background at the 
near-infrared wavelength bands. 

This work is based on observations with AKARI, a JAXA project with the participation 
of ESA. We thank the AKARI IRC team for their encouragement and helpful discussions. 
H.M.L. and M.G.L. were supported by NRF grant No. 2012R1A4A1028713. This research 
made use of Montage, funded by the National Aeronautics and Space Administration’s 
Earth Science Technology Ofhce, Computation Technologies Project, under Cooperative 
Agreement Number NCC5-626 between NASA and the California Institute of Technology. 
Montage is maintained by the NASA/IPAC Infrared Science Archive. 


Facility: Akari (IRC). 
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Fig. 1.— Mapping strategy of N2 band in the NEP deep field, which consists of central 
part, inner annulus, and outer annulus. Solid line box represents FOV of each pointed 
observation that covers about 10' x 10'. Solid line boxes with thick borders in the outer 
annulus are the area used in this study, whereas a dotted line c ircle wit h a diameter of 


10' is the area used in the previous study, i.e., the Monitor field fjPaper ll ). A dotted line 


box shows the position where the simulation for the clustering of unresolved faint galaxies 
is performed (see Section | 


by using the CFHT WIRCam Kg band image, which covers 
38' X 38'. Background image is the Improved Reprocessing of the IRAS Survey (IRIS; 
Miville-Deschenes fc Lagachel 120051 ) 25 pm map whose FOV corresponds to 64' x 63', which 
is retrieved from http://irsa.ipac.caltech.edu/data/IRIS/. Bright (white) region represents 


a planetary nebula NGC 6543. 


















37 



65 70 75 80 85 Mean of row (nW nn’^ sr"') 65 70 75 80 85 


Fig. 2.— Example images (pointing ID; 2110163) of N3 band are introduced to explain the 
correction for MUXbleed problem. Left image and graph describe a state before MUXbleed 
correction, whereas right parts describe a state after correction. Open circle represents the 
average background value of each row. Dotted line represents the position of MUXbleed, 
and we assume that the effect of MUXbleed reaches to the position where dot-dashed line 
is presented. Dashed line in the left graph is a htted line of data points that are regarded 
not to have been affected by MUXbleed, i.e., it is the base line for the correction. Solid 
line in the left graph is a htted line of data points that are thought to have been affected 
by MUXbleed. Correction for MUXbleed problem was performed by using the difference 
between solid line and dashed line in the left graph. 
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Fig. 3.— Sky brightness of each image for two bands. Open squares and open circles 
correspond to 2.4 and 3.2 pm, respectively. Solid lines are sinusoidal htting of the sky 
brightness as a function of heliocentric ecliptic longitude of Earth. Error bar represents Icr 


error. 
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Fig. 4.— Mosaic images of the NEP deep field for two bands in the equatorial coordinates. 
Top: hnal states of the mosaic images. White color represents masked regions. The small area 
denoted by dashed line box is the region where the effect of MUXbleed is small compared 
to other regions in the mosaic image regardless of the correction for MUXbleed problem. 
Bottom: smoothed mosaic images whose pixel values are replaced by non-weighted averaging 
of neighbor pixels within a circular shape hlter with a diameter of 100" to show large scale 
structures. The size of the mosaic image (solid line box) is 36' x 25' and 35' x 23' for 2.4 
and 3.2 //m, respectively. Note that structures in the mosaic images include artihcial as well 
as real ones in the sky. Mean pixel value of each mosaic image is adjusted to zero. 
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Fig. 5.— Top: fluctuation spectra, [q^P 2 {q)/2 'kY/‘^^ of the mosaic image (filled circles) and 
subset analysis (open circles), respectively. Dotted line represents the expected slope from 
random noise. Bottom; fluctuation spectra of the mosaic images after quadratically subtract¬ 
ing those of subset analysis. P 2 {(l) and q represent power spectrum and angular wavenumber, 
respectively. Error bars show la error. 
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Fig. 6.— Fluctuation spectra of the mosaic images before (asterisks) and after (open circles) 
MUXbleed correction for 2.4 and 3.2 pm, respectively. Solid line is fluctuation spectrum 
at the small area where the effect of MUXbleed is negligible compared to other regions in 
the mosaic image. We can see that the fluctuation spectrum after the MUXbleed correction 
is very similar to that of the area where MUXbleed is negligible up to ~300". Error bar 
represents la error. 
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Fig. 7.— Fluctuation spectra of the subset analysis with alternating pointed observations 
(open triangles) and alternating images (asterisks) for 2.4 and 3.2 pm, respectively. Error 
bar is Icr error. Dotted line represents the expected slope from random noise. 
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Fig. 8.— Here we show 13 images (N2 band) with pointing IDs (numbers above images) 
that are produced by subtracting the stacked image of each pointed observation of the 
Monitor held separately from the entire stacked image of the Monitor held. The last image 
is the stacked image of 13 images with pointing IDs. Images are smoothed by non-weighted 
averaging of neighbor pixels within a circular shape hlter with a diameter of 50" to emphasize 
structures. Black color represents masked regions. The diameter of each image is about 
550". Images are aligned with the equatorial coordinates. Mean pixel value of each image is 
adjusted to zero. Scale range corresponds to ±5cr of 13 images with pointing IDs. 
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Fig. 9.— Fluctuation spectra of the individual images (13 images with pointing IDs in 
Figure IHl solid lines) and the stacked image (the last image in Figure El open circles) for 2.4 
and 3.2 /im, respectively. Filled circles represent mean of the solid lines. Dotted line and 
dashed line represent the expected slope from random noise for each case. All fluctuation 
spectra are estimated by using images before smoothing. Error bars show Icr error. Error 
bars of the solid lines are not plotted to avoid confusion, while those of the hlled circles are 
too small to be seen. 
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Fig. 10.— Test of the mode coupling matrix for 2.4 and 3.2 /mi, respectively. Squares are 
seed fluctuations for simulations, and triangles are fluctuation spectra of the simulations after 
adopting mask of the mosaic image. Solid lines are recovered fluctuations by the inverse of 
the mode coupling matrix, M~^[q). Symbols and error bars represent mean and standard 
deviation of 500 simulations. 
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Fig. 11.— The beam transfer function of the NEP deep held. Symbols represent power 
spectra of PSFs that are normalized to be 1 at large angular scale. Squares and triangles are 
results for 2.4 and 3.2 /mi, respectively. Error bars represent a standard deviation of PSFs 
in pointed observations. 
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Fig. 12.— The map-making transfer function by the mosaicking procedure. It shows a ratio 
of remade fluctuation spectra to original ones. Squares and triangles are results for 2.4 and 
3.2 /im, respectively. Symbols and error bars represent mean and standard deviation of 
1000 simulations. We used observed fluctuation spectra at angular scales larger than 10" for 
further analysis because suppression of the fluctuation becomes signihcant at angles smaller 
than 10". 
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Fig. 13.— Top: the final fiuctuation spectra (filled circles) for 2.4 and 3.2 pm, respectively. 
Dashed line represents the contribution of shot noise that is estimated by simulations (see 
text for details). Error bar of the filled circle represents quadratic sum of all error components 
that are presented in the bottom panel. Bottom: all components of the error are plotted 
separately by various lines: la error of Fourier elements in each bin (long dashed), fiat 
field error (dotted), residuals of detected objects (dashed), and errors related with the mode 
coupling matrix (triple dot dashed), the beam transfer function (thin solid), and the map¬ 
making transfer function (dot dashed). Thick solid line is quadratic sum of them. 
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Fig. 14.— Left: correlations between the near-infrared mosaic images and the far-infrared 
mosaic image at 90 /im taken by AKARI Error bars represent the standard deviations 
of the noise images. Linear htted lines are shown as solid lines. The slopes of htted lines 
are almost zero (10“^ or so), indicating that there are no correlations between near- and 
far-infrared bands. Right: flnctuation spectra of the far-infrared mosaic image at both the 
Monitor held (asterisk) and the NEP deep held (open triangle). Because shot noise is a 
dominant component at small angular scale, we focused on huctuations at angles larger than 
100". Error bars show la error. 
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Fig. 15.— Simulation to estimate the clustering of faint objects. Filled circle is fluctuation 
spectrum of the mosaic image of 2.4 pm, and dashed line is shot noise of 2.4 pm. Solid 
line shows the contribution of the clustering of faint objects, which is estimated by using 
sources fainter than 21.5 AB mag at Kg band of CFHT WIRCam. Mean contribution of the 
clustering of faint objects, i.e., subtraction of the dashed line from the solid line in quadrature, 
at angles larger than 300" is 0.059 nW m“^ sr“^. Error bar of the solid line represents Icr 
error, and that of the dashed line means a standard deviation of 100 simulations. 
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Fig. 16.— Mean excess fluctuations at angles larger than 100". Open squa re repres ents the 


result of the NEP deep held, an d open circle is that o 


triangle is the result of Spitzer ( 


hlled diamond is that of GIBER (jZemcov et ah 


K^ashlinskv et al 


the Monitor held fjPaper ll ). Open 


20051) at a few hundreds arcseconds, and 


2014J ). Solid line represents a Rayleigh-Jeans 


spectrum (~A ^), which is htted by the open circle. 
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Fig. 17.— Excess fluctuation of the Monitor field (filled circles) and the NEP deep held 
(open squares). Their error bars are estimated by the error propagation formula with the 
error of the shot noise. Solid lines represent observational result of Spitzer data at 3.6 pm 


flKashlinskv et al. 


20121) . which is adjusted by using Rayleigh-Jeans spectrum (~A ^) for 


comparison with the result of AKARI. 
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Table 1. Sinusoidal fitting of zodiacal light 


Band 

a 

b 

c 


(nW sr“^) 

(deg.) 

(nW sr“^) 

N2 

11.70 ± 0.010 

66.20 ± 0.089 

113.92 ± 0.016 

N3 

6.70 ± 0.004 

61.62 ± 0.055 

71.26 ± 0.005 
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Table 2. Basic parameters of the mosaic image of the NEP deep held 


Item 

2.4 pm 

3.2 pm 

R.A.(J2000) DEC.(J2000) 

17'^55”^24" 

66°37'32" 

Number of pointed observations 

17 

17 

Number of images 

75 

68 

Mean observation time (minutes) 

12.0 

10.9 

Mean number of stacked images 

11 

10 

Limiting magnitude (AB mag) 

21.4 

21.8 

Percentage of remaining pixels (%) 

39.6 

39.5 

Average sky brightness (nW m“^ sr“^) 

113.92 ± 0.016 

71.26 ± 0.005 






